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A Method  for  Determining  the  Point  of  Llfi-Off 
and  Modified  Trajectory  of  a Ground-Based 
Heated  Turbulent  Planar  Jet  in  a 
Co-Flowing  Wind 


l.  IMROmCTION 

As  part  of  an  investigation  to  aid  in  the  development  of  an  operational  Warm 
Fog  Dispersal  System  (WFDS),  experimental  and  theoretical  studies  have  been 
made  of  the  characteristics  of  ground  based  heated  jets  with  varying  combinations 
of  heat  and  thrust  under  different  wind  conditions.  The  numerical  results  obtained 
from  these  investigations  will  be  utilized  to  help  develop  the  control  system  for  the 
WFDS.  As  part  of  the  theoretical  study,  Klein  and  Kunkel^’  described  methods 
for  calculating  the  characteristics  of  planar  and  round  jets  in  a co-flowing  wind; 
that  is,  in  the  same  direction  as  the  jet.  The  characteristics  of  heated  jets  in  a 
counter-flowing  air  stream  will  be  present  in  subsequent  reports. 

In  the  development  of  the  previous  models,  the  effect  of  the  ground  was  neg- 
lected and  the  buoyant  thrust  of  the  air  upon  the  jet  was  assumed  to  start  acting  as 
soon  as  the  heated  gas  exited  from  the  nozzle.  The  calculated  trajectories  conse- 
quently start  to  rise  above  the  ground  a short  distance  from  the  jet  nozzle.  Field 


(Received  for  publication  1 February  1977) 

1.  Klein,  M.M.,  and  Kunkel,  B.A.  (1975)  Interaction  of  a Buoyant  Turbulent 
Planar  Jet  With  a Co-Flowing  Wind,  A FCRL-TR -75-0368. 


results  described  by  Kunkel^  show,  however,  that  the  jet  trajectories  generally 
hug  the  ground  for  a considerable  distance  before  the  effect  due  to  buoyancy  is 
exhibited.  This  ground  effect  appears  to  stem  from  the  development  of  a regime 
of  recirculating  flow  in  which  the  pressure  is  less  than  that  of  the  surrounding  air, 
Conway.  ^ A complete  description  of  the  original  model  and  its  derivation  for  the 
planar  case  is  contained  in  Klein  and  Kunkel.  ^ The  present  report  extends  this 
model  to  take  account  of  the  ground  effect. 

Although  desirable  to  be  able  to  calculate  the  point  of  lift-off  from  first  prin- 
ciples, it  was  felt  that  such  an  approach  would  be  lengthy  and  complex.  Further, 
in  view  of  the  incomplete  knowledge  of  the  details  of  the  flow,  the  solutions 
obtained  might  not  be  accurate  enough  for  utilization  in  the  design  of  a control 
system  for  the  WFDS.  It  was,  therefore,  decided  to  analyze  the  trajectory  data 


available  from  the  field  experiments  to  see  if  the  point  of  lift-off  could  be  deter- 

3 

mined  from  the  physical  parameters  defining  the  flow.  Since  Kunkel  found  that 


the  jets  merge  at  about  2.  5 times  the  combustor  soacing,  that  is,  a short  distance 
down  stream,  the  analysis  was  made  for  the  planar  jet  case. 

An  examination  of  the  data  shows  that,  as  ant  cipated,  the  point  of  lift-off  is 
delayed  by  increasing  the  velocity  of  the  jet  and  hastened  by  increasing  the  jet 
temperature.  An  obvious  possible  parameter  would  be  the  ambient  Froude  number, 
that  is,  the  ratio  of  inertial  to  convective  forces  in  the  initial  jet.  However,  the 
lift-off  point  must  also  depend  upon  windspeed  and  should  be  intimately  related  to 
the  temperature  and  velocity  at  this  location.  The  analysis  was,  therefore,  made 
in  terms  of  the  local  Froude  number  at  the  point  of  lift-off. 


2.  \wnsi? 


The  point  of  lift-off,  x^,  axial  velocity,  u 
over  ambient,  AT 


and  axial  jet  temperature  excess 


c'  •"  "m’ 

are  indicated  schematically  in  Figure  1,  The  local  Froude 


number,  Fj,  at  the  point  of  lift-off,  may  be  written  as. 


u^  T 
m a 


1 


(1) 


Kunkel,  B.A.  (1975)  Heat  and  Thrust  Requirements  of  a Thermal  Fog 
Dispersal  System,  AFCRL-TR-75-04V2. 


4.  Conway,  Arthur  (Editor)  (1975)  A Guide  to  Fluidics,  MacDonald  & Co. 
(Publishers,  Ltd),  London,  pp.  2-7, 
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> 

AMBIENT  WIND.Ug 


Figure  1.  Schematic  Representation  of  Jet  and  Lift-Off  Point 


where  T is  the  ambient  temperature,  g the  gravitational  constant,  and  the 

initial  jet  width.  The  functional  relation  f,  if  it  exists,  between  and  may 

then  be  expressed  as 


Xc  = f(F,) 


(2) 


It  is  convenient  to  write  F^  in  the  form  (the  subscript  zero  refers  to  initial  jet 
values) 


F,  = 


2 2 

o m 


1 ,2 

O \1  « Vk 

o g b 


at: 


(3) 


o AT, 


;_1 


and  express  AT^  in  terms  of  the  velocity  incr  lent  z by  (see  Kunkel  and  Klein 
for  more  detail) 


AT 


= kz 


(4) 


where 


u - u 
m a 


u - u 
o a 


(5) 
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and  u is  the  wind  velocity,  and  k a coefficient  which  varies  slowly  with  jet  and 
wind  speeds. 

The  Froude  number  may  now  be  written  as 


(6) 

■) 

where  m is  the  ratio  of  wind  velocity  to  initial  jet  velocity. 

For  purposes  of  plotting,  it  is  convenient  to  introduce  a new  number  in 
which  the  constant  and  nearly  constant  terms  have  been  omitted,  that  is. 


(7) 


To  see  if  Gj^  is  a proper  parameter  for  correlation  of  the  data,  a plot  was  made  of 
X against  Gj  and  a smooth  curve,  estimated  by  eye,  drawn  through  the  points 
(Figure  2).  An  examination  of  Figure  2 shows  that,  while  most  of  the  points  lie 
close  to  the  curve,  several  of  the  points  - identified  by  tick  marks  - lie  distinctly 
below.  Since  these  points  occur  at  the  higher  values  of  m(m  ' 0.  1),  it  appears 
that,  as  the  windspeed  increases,  the  point  of  lift-off  decreases.  This  behavior, 
which  is  opposite  to  the  expected  increase  of  x^  with  windspeed,  may  be  explained 
in  terms  of  the  regime  of  low  pressure  recirculating  air  described  by  Conway. 

The  decrease  in  pressure  of  the  recirculating  air  will  be  determined,  to  a large 
extent,  by  the  increase  in  velocity  of  the  air  mixing  with  the  jet.  As  the  windspeed 
increases,  a smaller  velocity  increase  is  obtained,  resulting  in  a weaker  pressure 
decrease  and  earlier  lift-off. 

To  take  into  account  this  behavior,  the  Froude  number  F^  in  Eq.  (1)  was 
modified  to 


F = 


u (u  - u,^)T„ 
m m a o 

gb„  AT 
“ o o 


resulting  in,  as  the  modified  form  of  Eq.  (6) 


F = 


(8) 


(9) 
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Figure  2,  Plot  ot  l.ift-Off  Point  Against  Parameter  Gj  for  Co-Flowing 
Test  Data.  Points  with  tick  marks  represent  higher  values  of  m(m-'0.  1) 
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Figure  3.  Plot  of  Lift-Off  Point  x^.  Against  Parameter  G for  Co-Flowing 
Test  Data.  Solid  curve  is  power  law  fit 


As  in  the  case  of  Fj^,  we  utilize  a number  G where 


G = 


% 

b AT  u 


o o o 


(10) 


in  which  the  constant  and  nearly  constant  terms  have  been  omitted.  A plot  of 

against  G,  Figure  3,  now  shows  that  all  the  points,  including  those  associated 

with  the  higher  values  of  m form  a reasonably  defined  curve.  It  may  be  noted  that 

if  we  had  modified  Eq.  (1)  with  the  term  u z rather  than  (u  - u ),  Eq.  (10)  would 

o m a 

lack  the  factor  (1  - m).  The  effect  of  this  term  is  almost  negligible  except  at  the 
higher  values  of  m,  where  it  results  in  a slight  improvement  of  the  data.  Addi- 
tional reliable  data  at  moderate  to  higher  values  of  m,  would  be  required  to  deter- 
mine more  precisely  the  dependence  of  the  lift-off  point  upon  m. 

The  data  points  in  Figure  3 can  be  well  fitted  by  the  power  law  curve  (for 
purposes  of  correlation,  distances  are  dimensionless  with  b^  the  reference  length) 


X 


c 


190.  0 


„0. 744 

u 


(11) 


The  velocity  u^  is  a function  of  x and  may  be  calculated  from  z in  Eq.  (5)  by 


u 

-jj2-  = z(  1 - m)  + m 


where,  at  x = x^,  z is  given  by  (see  Klein  and  Kunkel^) 


2 ^ 1 J_ 

^ " T 1 - m A„  ' 

Si  A 


(12) 


(13) 


(14) 


(15) 


Aj  and  A2  are  profile  constants  and  x^  is  the  initial  jet  position.  Although  Eq.  (11) 
can  be  solved  numerically  for  x^,  a formal  solution  is  not  available  because  of  the 
nonintegral  exponent  occurring  in  G.  However,  an  explicit  solution  can  be  obtained 
if  the  variables  x^  and  G are  changed  to  q and  p,  where 


12 


(16) 


and 

u^(l  - m) 
„ o 
P AT  b 
o o 


i 

I 

(17) 


Note  that  the  parameter  p is  independent  of  and  can  be  immediately  calculated 
for  a given  set  of  flow  conditions.  Since  the  exponent  of  G in  Eq,  (11)  does  not 
differ  greatly  from  unity,  we  may  expect  good  correlation  in  a plot  of  q against  p. 
This  anticipated  result  is  confirmed  in  Figure  4 where  the  scatter  of  the  plotted 
points  is  seen  to  be  small.  A good  fit  to  the  data  is  provided  by  the  power  law 
curve 

q = 253.  64  p°-®°  (18) 


with  a correlation  coefficient  of  0.  98. 


Figure  4.  Plot  of  Parameter  q Against  Parameter  p for  Co-Flowing  Test  Data. 
Solid  curve  is  power  law  fit 
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To  obtain  the  lift-off  point  for  a given  value  of  the  parameter  q,  we  utilize 
Eqs,  (12)  Eind  (13)  to  express  q in  Eq.  (16)  as  a function  of  x^,  3rielding 


X V/2 


(19) 


where 


Cj  = (1  - m)  a 


and 


C2  = o - m 


Rationalization  of  Eq.  (19)  leads  to  the  cubic 


3 2 

X + a„x  + a,  x„  + a = 0 
c 2 c 1 c o 


(20) 


where 


^2  = 2c2  q 

®1  * ^*^2  "“2^1^  " -m(2<yCj  - m)  q^ 

2 2 

A study  of  Eq.  (20)  for  the  range  of  parameters  encountered  in  the  field  data  yields 
two  complex  roots  and  one  real  root.  The  real  root  is  given  by  Abramovitz  and 
Stegun^  (p  17), 

Xc  = Si  +S2  - 82/3  , (21) 


5.  Abramovitz,  M. , and  Stegim,  I.  (Editors),  (1964)  Handbook  of  Mathematical 
Functions,  U.  S.  Government  Printing  Office,  Washington,  D.  C. 
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where 


1/3 


[r  + «’  ♦ 

2 p 2 T 8 3 

= 2m(m  - oCj)(2aCj  - m)q  + 3x^Cj  I + - ac^) 

4[< 


I 4o^  c j - m (2a  Cj^  - m)j  . 


3 

q 


Since  the  jet  characteristics  at  a given  location  are  independent  o/  the  vertical 
position,  y,  the  method  of  calculation  presented  in  Klein  and  Kunkel''  requires  only 
minor  modifications  to  take  account  of  the  ground  effect.  The  buoyant  force  is 
taken  as  zero  up  to  the  lift-off  point  x^,  at  which  point  it  is  assumed  to  start  acting. 
The  jet  characteristics  at  this  point  are  then  calculated  and  utilized  as  the  initial 
conditions  for  the  subsequent  buoyant  motion  of  the  jet.  Of  interest  for  fog  dis- 
persal work,  in  addition  to  the  jet  trajectory,  is  the  temperature  distribution  along 
the  jet,  which  is  easily  obtained  from  Eqs.  (4)  and  (13). 


3.  RESULTS  AM)  OISaiSSIO\ 

The  calculated  trajectories  with  and  without  ground  effect  and  the  correspond- 
ing experimental  curves  are  presented  in  Figure  5.  These  results  cover  most  of 
the  range  of  initial  jet  speed,  initial  jet  temperature,  and  ambient  wind  encountered 
in  the  field  tests.  The  calculated  trajectories  are  in  fair-to-good  agreement  with 
the  experimental  curves,  generally  tending  to  be  somewhat  higher,  especially  at 
very  low  windspeeds.  In  addition,  several  of  the  calculated  lift-off  points  are  not 
in  good  agreement  with  the  experimental  values,  a result  due  to  the  fact  that  several 
of  the  experimental  lift-off  points  did  not  lie  close  to  the  correlation  curve  of  q vs. 
p (Figure  4).  This  behavior  in  the  lift-off  point  may  be  due  to  an  instability  of  the 
air  near  the  ground  with  respect  to  fluctuations  in  the  ambient  wind,  under  certain 
flow  conditions.  In  test  16-9,  the  experimental  curve  shows  a flattening  tendency 
at  large  jet  distances  not  exhibited  by  the  calculated  trajectory.  An  examination 
of  the  upper  level  wind  read  by  the  wind  sensor,  shows  that  it  is  almost  double  the 
ambient  wind.  A correction  to  the  later  portion  of  the  jet  trajectory,  utilizing  the 
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Figure  5,  Comparison  of  Calculated  Trajectories,  for  Several  Tests,  With  Corre 
spending  Experimental  Curves 
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Figure  5.  Comparison  of  Calculated  Trajectories,  for  Several  Tests,  With  Corre- 
sponding Experimental  Curves  (Cont) 


proper  jet  influenced  upper  level  wind,  brings  it  into  good  agreement  with  the 
experimental  curve.  This  flattening  out  of  the  trajectory  is  exhibited  principally 
by  those  tests  in  which  the  ambient  wind  was  strong  enough  to  keep  the  trajectory 
from  becoming  very  steep. 

Typical  calculated  and  experimental  temperature  distributions  along  the  axis 
are  presented  in  Figure  6 for  two  initial  jet  temperatures.  In  general,  the  calcu- 
lated curves  decrease  at  a somewhat  slower  rate  than  those  obtained  experimen- 
tally. This  results  in  temperatures  that  are  too  high  at  large  jet  distances,  with 
the  effect  being  greater  for  the  case  of  higher  initial  jet  temperature.  Since  the 
jet  is  strongly  inclined  to  the  vertical  at  large  jet  distances,  the  value  of  the  mix- 
ing coefficient  c — valid  for  a co-flowing  nonbuoyant  jet  — may  be  inaccurate  in  this 
region.  In  addition,  the  physical  trajectory  will  be  deformed  by  the  dynamic  pres- 
sure exerted  by  the  ambient  wind  upon  the  jet  when  the  angle  between  the  two  flows 
is  large.  A more  accurate  model,  incorporating  these  and  other  possible  mod- 
ifications, is  beyond  the  scope  of  the  present  investigation. 


6.  Abramovich,  G.  N.  (1963)  The  Theory  of  Turbulent  Jets,  The  MIT  Press, 
Cambridge,  Mass.,  pp.  541-556. 
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4.  PROCEDI  RE  FOR  RAPID  DETERMINATION  OF  LIFT-OFF  POINT 


For  use  in  field  work,  it  is  desirable  to  be  able  to  rapidly  estimate  the  lift- 
off point  of  the  jet.  An  examination  of  Section  2 shows  that  depends  explicitly 
upon  m and  and  implicitly  upon  b^  and  u^  through  the  parameter  p.  The 
solution  for  x may,  therefore,  be  presented  in  fairly  compact  form  without 
explicit  reference  to  b^  and  u^  provided  p is  used  as  the  independent  variable. 
Calculations  of  x^  have  accordingly  been  made  for  a range  of  input  data,  and  plotted 
in  Figure  7 against  p,  with  m and  as  parameters.  As  seen  from  the  figure, 
the  lift-off  point  is  much  more  sensitive  to  m than  to  ^T^,  the  weak  dependence 
upon  the  latter  becoming  very  small  when  m is  increased  to  0.  2. 

For  a given  set  of  flow  conditions,  the  value  of  x^  is  easily  obtained  from 
Figure  7.  As  sui  illustrative  example  we  take: 

u = 16  m/sec  , u = 0. 9 m/sec 

T = 175  K°  , b = 0, 01  m . 

o * o 

The  wind  parameter  m is  given  by 

m = ^=  0.  056  . 

u lb 
o 

The  parameter  p may  now  be  calculated  from 

u^  (1  - m) 

_ _ _o 

P ■ b 

o o 


to  yield 

256(0.944) 

P ■ 0.01(175) 

= 138. 1 . 

From  Figure  7,  by  interpolation 

X = 1520,  in  units  of  b 
c ’ o 

= 15.  2 m . 
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5.  SI'MMARV  AM)  CONCLUSIONS 


The  method  of  calculating  the  dynamic  characteristics  of  a ground  based 
turbulent  heated  planar  jet  has  been  extended  to  take  account  of  the  tendency  of  the 
jet  to  remain  attached  to  the  ground  for  a considerable  part  of  the  jet  trajectory. 
The  lift-off  points,  determined  from  the  experimental  data,  have  been  analyzed 
and  correlated  in  terms  of  the  local  Froude  number  at  the  lift-off  point.  From 
this  analysis,  a procedure  has  been  developed  for  determining  the  lift-off  point 
when  the  flow  conditions  for  the  jet  — that  is,  ambient  wind,  initial  jet  velocity,  and 
temperature  - are  given.  The  new  jet  trajectory  may  now  be  calculated,  with 
only  a minor  modification  of  the  original  procedure.  The  calculated  trajectories, 
and  temperature  distributions  along  the  jet  are  in  fair-to-good  agreement  with  the 
corresponding  results  from  the  field  data.  The  present  method  of  calculating  the 
jet  trajectories  appears  accurate  enough  to  be  utilized  in  the  development  of  a con- 
trol model  for  the  fog  dispersal  system. 
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